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ABSTRACT 



The majority of the heavy elements produced by stars 2 billion years after the Big 
Bang (redshift z ~ 3) are presently undetected at those epochs. We propose a solution to 
this cosmic 'missing metals' problem in which such elements are stored in gaseous halos 
produced by supernova explosions around star-forming galaxies. By using data from the 
ESO/VLT Large Program, we find that: (i) only 5%-9% of the produced metals reside 
in the cold phase, the rest being found in the hot (T = ifj 5 - 8-6 - 4 K) phase; (ii) l%-6% 
(3%-30%) of the observed C IV (O VI ) is in the hot phase. We conclude that at z }t 3 
more than 90% of the metals produced during the star forming history can be placed in 
a hot phase of the IGM, without violating any observational constraint. The observed 
galaxy mass-metallicity relation, and the intergalactic medium and intracluster medium 
metallicity evolution are also naturally explained by this hypothesis. 

Subject headings: (galaxies:) intergalactic medium - supernovae: general - galaxies: 
stellar content - stars:early-type 



1. Motivation 

In its original formulation (Pettini 1999), the 'missing metals' problem was stated as follows. 
Studies of the comoving luminosity density of distant galaxies allow us to trace the cosmic star 
formation density (or history, SFH), p*(z), up to redshift z max ~ 7. Assuming an initial mass func- 
tion of such stars (IMF), one can compute the specific fraction of heavy elements ('metals') they 
produce, y, and derive the metal production rate pz(z) = yp+(z), whose integral from z max gives 
the density of cosmic metals in units of the critical density, 0, s / h , at any given z. Early searches in 
cosmic structures for which the metal/baryon mass ratio 1 (metallicity, Z = Qz/^b) can be derived 
either via intergalactic gas quasar absorption line experiments (Damped Lya Absorbers [DLAs] 
or the Lya 'forest') or through direct spectroscopic studies (Lyman Break Galaxies [LBGs]) have 
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found that only Q°^ s £ 0.20Q s J h is stored in these components, i.e., the large majority of the metals 
are 'missing'. An analogous missing metal problem is also found by considering in a self-consistent 
manner the star formation rates and metallicities of DLAs alone (Wolfe et al. 2003, Prochaska et 
al 2003). Newly available high-quality data allow a more precise analysis of the problem. 



2. Stating the problem 

To re-evaluate £l s J h we use the most recent SFH compilation (Bouwens et al. 2004) corrected 
upwards for the effects of dust obscuration by the prescribed (Reddy & Steidel 2004) value of 4.5 
at z > 3 and adopt y = 1/42. Integration of p z (z) to z = 2.3 yields Sl s / h = 1.84 ± 0.34 x 10~ 5 . 
Where should we look for these metals ? 

The most obvious locations are the galaxies used to derive the SFH, i.e., LBGs. These are char- 
acterized (Reddy & Steidel 2004) by a mean comoving number density of 6 x 10~ 3 h 3 Mpc~ 3 and 
Z = 0.6Z Q = 0.0113. Stellar masses can be constrained only by assuming a range of star formation 
histories of the form SFR(t) oc exp(— t/r) and therefore they are somewhat uncertain. According 
to Shapley et al. 2005, they should be in the range 0.6 — 6 x 10 10 Mq. Assuming the best fit value 
= 2 x 1O 1O M , we get d h z 9 = 3.4 x lO^M^io « 0.18^ //l - If metals are not in LBG stars or gas, 
they could be in DLAs or the IGM. The metal content of DLAs is derived by noting that (Rao & 
Turnsheck 2000, Prochaska & Wolfe 2000) at z ss 3 their neutral (« total) gas density Qf a = 10~ 3 
and metallicity Z = 3.8 x 10~ 4 combine to give ftf° = 3.8 x 10~ 7 < Uf s < Cl s J h ; we can therefore 
neglect their contribution. In the following, we correct for LBG contribution by (re)defining the 
cosmic density of missing metals £l s J h = Q s J h — S7^ 9 ' . 

Hence, the missing metals should be found outside the main body of galaxies or DLAs. There are 
essentially two possibilities: (a) they could reside in the Lya forest, or (ii) in the halos of the galax- 
ies that have produced and ejected them. Note that the distinction between these two components 
is somewhat ambiguous. Our working definition is that galactic halos are gravitationally bound 
structures around galaxies; they are special sites affected by galactic winds. 

The most widely studied tracers of Lya forest metallicity are C IV and O vi absorption lines. 
The fraction of the critical density p c contributed by a given element ion, Ei, of mass the residing 
in the Lya forest is given by 

where AA(z + ,z_) = f* + dz(l + z) 2 E(z)~ l is the absorption distance (Bahcall & Peebles 1969), 
with E(z) = [J7a + ^m(l + z) 3 ] 1 / 2 ; sums are performed over all the redshift intervals z_ < z < z + 
over which an ion column density Ne { is detected. To determine ^qvi anc ^ ^CIV' we use ^ a * a f rom 
the ESO VLT Large Program 2 (Bergeron et al. 2002) which provides high S/N, high resolution 
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spectra of an homogeneous sample of 19 QSOs in 1.7 < z < 3.8; ^ovi * s curren tly available for 
four LP lines of sight (Bergeron et al. 2002, Bergeron Sz Herbert-Fort 2005), for which we find 
S7qyj = 1.3 x 10~ 7 ; two other recent estimates (Simcoe, Sargent & Rauch 2004, Carswell, Schaye 
&; Kim 2002), give ^ovi = (1-1) 0-9) x 10~ 7 . We adopt the sightline-weighted mean of the three 
values allowing for the largest error, ^qvi = ^ ^ ^ x F rom the C IV absorber distribution 

(Aracil et al. 2004, Scannapieco et al. 2005) in the column density range 12 < logAciv < 16 
we obtain = 7.5 ± 2.2 x 10~ 8 (statistical error). This value is about two times higher than 

previous determinations (Songaila 2001; Simcoe, Sargent & Rauch 2004, Schaye et al 2003) which 
could not account for the contribution of strong (log C IV > 14) absorption systems. Combining the 
average measured Nqw-Nhi and Noyi-Nhi correlations (Aracil et al. 2004) with the measured 
distribution of weak HI absorbers (Petitjean et al 1993), we have checked that systems with 
log iVciv < 12 contribute less than 1%, well within the quoted error. For a (meteoritic) solar carbon 
logarithmic abundance (in number) Ac = 8.52 with respect to hydrogen (Ah = 12), we conclude 
that only a fraction Q,Q^/Q s J h = 2.4 x 10~ 2 of the produced carbon is observed in the C IV state. 
Repeating the procedure for O (Ao = 8.83), gives a ratio f^yj/fi^f 1 = 1.3 x 10~ 2 . To account for 
all uncertainties above, we will consider values in the range i.4 x io- 2 < n l ^/n s J h < 4.0 x io- 2 
and 8.1 x IO" 3 < fig^/fi^ < 2.1 x IO" 2 . 

We now determine the physical conditions of the gas hiding the missing C and O. Numerical 
simulations (Dave et al. 2001) suggest that the intergalactic medium [IGM] might be a two-phase 
system made by a cool (T c ~ 10 4-4 ' 5 K), photoionized phase, and a hot, collisionally ionized one. 
We impose the following conditions separately for each ion (C IV ,0 VI ) and element (C,0): (1) the 
observed ionic abundance is the sum of the abundances in the two phases; (2) the SFH-produced 
element abundance is the sum of the element abundances in the two phases; (3) the elements 
are in the same abundance ratios in the two phases. More explicitly, these conditions can be 
mathematically expressed as 

fc^c + fc^c = ^civ (2) 
fo^o + fo^o = ^ovi (3) 

n c + n h c = % fh (4) 
n c o + ^o = ^o h ( 5 ) 

n c c -AQ c = (6) 
QJc-A^o = (7) 

After some simple algebra, the above equations reduce to: 

ci l v a io s f h — fh ct l v a /r> s f h fh 
Si civ/ 5i c* Jc _ iL ovi/ iL o J( 
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where // = j\ (Aj,Tj,Uj) is the ionization correction for the considered ion (C IV or O VI ) of a 
given element (i = C,0) in the cold or hot phase, (j = c,h), respectively; the overdensity, Aj and 
temperature, Tj, of the two phases are the unknowns of the problem; finally, A is the abundance 
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ratio of the two elements. We complement these conditions by further imposing that the pressure of 
the cool phase does not exceed that of the hot one and assuming a temperature-density relation for 
the cold phase T = TqAJ, (with To = 2 x 10 4 K and 7 = 0.3), as inferred from the Lya forest data. 
The value of the photoionization parameter, Uj = n^/nj, is fixed by the ionizing photon density 
n 7 of the assumed UV background spectrum (Haardt & Madau 1996) shifted so that the intensity 
at 1 Ryd is J u = 0.3 x 10~ 21 erg s" 1 Hz -1 = 0.3 J21, corresponding to a hydrogen photoionization 
rate T = 0.84 x 10~ 12 s _1 = 0.84 Ti2, in agreement with Bolton et al. 2005. Finally, we warn that 
deviations from solar abundances might be possible, and indeed there are hints that oxygen might 
be overabundant (Telfer et al. 2002; Bergeron et al. 2002); here we neglect this complication. 



3. A possible solution 

By solving eq. 8, we obtain the results plotted in Fig. 1. The hot phase is characterized by a 
wide density range, logA^ > 0.4 and a restricted temperature range 5.8 < logT^ < 6.4. We find 
that: (i) only 5%-9% of the produced metals reside in the cold phase, the rest being found in the 
hot phase; (ii) l%-6% (3%-30%) of the observed C IV (O VI ) is in the hot phase. We conclude 
that more than 90% of the metals produced during the star forming history can be placed in a 
hot phase of the IGM, without violating any observational constraint. To further constrain the 
hot phase parameter range, we have searched in the LP C iv line list for components with large 
Doppler parameters. We find no lines with bcw > 26.5 km s _1 , corresponding to logT^ > 5.7; this 
result seems to exclude the high density and temperature region of the allowed parameter space in 
the middle panel of Fig. 1. We checked that the above findings are insensitive to variations of Ti2 
of ±50%; however, O vi /C TV ratios in the cold phase might depend on the UVB shape around 4 
Ryd. 

The derived values of and A^ are suggestive of regions likely to be found around galaxies; more- 
over, 10 6 K gas temperature would have a scale height of > 10 kpc, hence it cannot be confined 
in the disk. To test this hypothesis we resort to cosmological simulations. As an illustration, Fig. 
2 shows the temperature and velocity structure in a 2D cut through the center of a simulated 
galaxy (we used the multiphase version [Marri & White 2003] of the GADGET2 code to simulate 
a comoving 10/i _1 Mpc 3 cosmic volume) at redshift z = 3.3; its total (dark + baryonic) mass is 
2 x 1O 11 M0, the star formation rate ~ 20Mq yr _1 . This galaxy has been selected to match LBG 
properties, but it is not unusual in the simulation volume. As often observed in LBGs, a strong 
galactic wind is visible, whose expansion is counteracted by energy losses due to cooling and grav- 
ity, and ram pressure exerted by the infalling gas. Infall is particularly effective at confining the 
wind into a broadly spherical region of physical radius ps 300 kpc, into which cold infalling streams 
of gas penetrate. Inside such wind-driven bubble the temperature (Fig. 2) is roughly constant 
T ps 10 6 K, whereas the density spans values of < log A < 5 [A(z = 3.3) = 1 corresponds to 
ps 2 x 10 -5 cm -3 ]. The cool phase is evident in the outer boundary of the bubble, where cooling 
interfaces arise from the interaction with infalling streams. Hence halos of LBGs seem to meet the 
requirements as repositories of missing metals. 
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Additional support for this conclusion comes from studies of the correlation properties of C IV and 
O VI absorbers (Pichon et al 2003, Aracil et al. 2004, Bergeron & Herbert-Fort 2005), which con- 
clude that: (i) O VI absorption in the lowest density gas is usually (about 2/3 of the times) located 
within ~ 300 — 400 km s _1 of strong H I absorption lines; (ii) the C IV correlation function is con- 
sistent with metals confined within bubbles of typical (comoving) radius « XAhr 1 Mpc in halos of 
mass M > 5 x 1O 11 M at z = 3. If each such objects hosts one bubble, the cosmic volume filling 
factor of metals is fz = 11%; it follows that halo metallicity is Q s J h / fz^b = O.165Z . A tempera- 
ture of log Th = 5.8 corresponds to H I (O vi ) Doppler parameters &hi = 102 (bovi = 25.5) km s _1 
and to Nqyi/Nui = 3; absorbers with logAovi = 13 are detectable for &ovi = 25.5 km s" 1 but 
the corresponding logA^ni = 12.4 ones for bm = 102 km s _1 are not. This raises the possibility of 
finding O VI absorbers without associated H I . 



4. Implications 

The scenario proposed leads to several interesting consequences. First, metals produced by 
LBGs do not seem to be able to escape from their halos, due to the confining mechanisms mentioned 
above. This is consistent with the prediction (Ferrara, Pettini & Shchekinov 2000) that galaxies 
of total mass A4 > 10 12 (1 + z)~ 3 / 2 M do not eject their metals into the IGM. Interestingly, 
the metallicity-mass relation recently derived from the SDSS (Tremonti et al. 2004) shows that 
galaxies with stellar masses above 3 x 10 10 M Q (their total mass corresponds to M for a star 
formation efficiency /* = 0.2) chemically evolve as "closed boxes," i.e., they retain their heavy 
elements. Second, the nearly constant (2 < z < 5, Z ~ 3.5 x 1O _4 Z ) metallicity of the low 
column density IGM (Songaila 2001) is naturally explained by the decreasing efficiency of metal 
loss from larger galaxies. Early pollution from low-mass galaxies allows a sufficient time for metals 
to cool after ejection; however, the majority of metals in LBG halos at lower redshifts are still 
too hot to be detected. Hence their contribution to the metallicity evolution of the IGM cannot 
be identified by absorption line experiments, which mostly sample the cool phase of the forest. 
Third, the rapid deceleration of the wind results either in a quasi-hydrostatic halo or in a 'galactic 
fountain' if radiative losses can cool the halo gas. In both cases this material is very poorly 
bound and likely to be stripped by ram pressure if, as it seems reasonable, the galaxy will be 
incorporated in the potential well of a larger object (galaxy group or cluster) corresponding to 
the next level of the hierarchical structure growth. Turbulence and hydrodynamic instabilities 
associated with this process are then likely to efficiently mix the metals into the surrounding gas 
within approximately a sound crossing time of ~ 1 Gyr, or Az ~ 0.5. If metals produced and 
stored in LBG halos by z = 2.3 end up in clusters, than the average metallicity of the intracluster 
medium is Zjcm = ^z ' I^icm = O.31Z , having assumed (Fukugita, Hogan & Peebles 1998) 
fi/CM = 0.0026/iyQ 1 ' 5 . Not only is this number tantalizingly close to the observed value at z = 1.2 
(Tozzi et al. 2003), but we also predict that little evolution will be found in the ICM metallicity 
up to z ss 2 as essentially all the metals that could have escaped galaxies during cosmic evolution 
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had already done so by this epoch. 
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Fig. 1. — Allowed hot phase gas temperature, Th, and overdensity, Ah, regions. Also shown are the corresponding 
isocontours of the Q.%/S~l s / h {upper panel), ilcw /^crv ("middle panel) and f2o VI /f2Qyj (bottom panel). 
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Fig. 2. — Bottom panel: Temperature map (physical Mpc units) of a plane through the center of a simulated 
starburst galaxy at z = 3.3 with properties typical of Lyman Break Galaxies; longest velocity vectors correspond to 
v — 150 km s _1 . Upper panel: ID cuts parallel to the horizontal axis and passing through the center of the map in 
the bottom panel showing the gas density (blue) and temperature (red) 



